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Abstract 

In light of recent measurements of the nucleon structure function in the small-x deep-inelastic regime 
at HERA and the consequently improved theoretical understanding of the quark distributions in this range 
of parton fractional momentum, we present new results for the neutrino-nucleon cross section at ultra high 
energies, up to £^ ~ 10^^ eV. The results are relevant to deep underground muon detectors looking for 
such neutrinos. For ~ 10'^" eV neutrinos, our cross sections are about a factor of 4 to 10 times larger 
than the previously reported results by Reno, Quigg and Walker. We discuss implications of this new 
neutrino-nucleon cross section for a variety of current and future neutrino detectors. 
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As the number of extant and planned experimental facilities for extra-terrestrial neutrino detection testifies, 
their importance in extending the particle physics frontier beyond the standard model is widely recognized 
[Q. At the highest energies, neutrinos are decay products of pions produced in cosmic accelerators and 
thus provide a direct window to the most energetic processes in the universe. Neutrino astronomy also has 
important advantages over gamma-ray astronomy. In the energy range 10^^ to 10^" eV, the expected air- 
cascade signature of the gamma-ray flux is dwarfed by the cosmic-ray background. In the same energy range 
the flux of neutrino by-products of the cosmic-ray interactions is very small, and the expected neutrino flux 
from cosmic accelerators should dominate. In addition, neutrino telescopes span a significant fraction of the 
sky at all times, in contrast to gamma-ray detectors which provide small angular coverage for a limited amount 
of time. 

Recent observations |Q of TeV gamma-rays from Mkn 421 and Mkn 501 have revived interest in studying 
mechanisms for producing high energy photons in Active Galactic Nuclei ( AGN) . If the observed photons are 
decay products of tt'^s produced in hadronic interactions in the disk surrounding the AGN, then AGNs are 
also powerful sources of ultra high energy (UHE) neutrinos . 

In general, fluxes of UHE neutrinos fall steeply with increasing energy, making the direct detection of these 
particles difficult. Cerenkov detectors, however, are capable of detecting upward-moving muons produced by 
the charged-current (CC) interactions of energetic neutrinos in the rock below the detector. The effective 
volume of the detector is thus enhanced in proportion to the range of the produced muon (typically several 
kilometers), which then may traverse the fiducial volume of the detector or stop therein. The expected 
event rate for these detectors is presently subject to two sources of uncertainty: a) the uncertainty in the 
knowledge of the incident neutrino flux, and, b) the calculation of the UHE charged current cross-section for 
1/^^ + N fj. + X, where is a nucleon. The latter stems from necessary extrapolations of the nucleon quark 
structure functions for very low parton fractional momentum x and large momentum transfer Q^. 

Using recent results from the ep collider at HERA Q on small- a; parton distributions, we have performed 
calculation of the neutrino-nucleon cross sections for neutrino energies up to 10^^ eV [|]. For 10^° eV neutrinos 
we have found the cross sections to be factor of four to ten times larger than previous estimates. Since UHE 
neutrinos are detected via neutrino- induced muons, our results translate into signiflcantly larger probability 
that neutrino-induced muon created in the rock surrounding the detector would reach the detector and thus 
larger downward muon rates. For upward-moving muons, the neutrino attenuation in the Earth depends on the 
charged-current and neutral-current cross sections. For larger neutrino-nucleon cross sections, the attenuation 
effect is stronger. Since it is the product of the muon probability and the neutrino attenuation that goes into 
the calculation of the upward-moving muon event rate, one has to incorporate these two competing effects, 
as elaborated below. 

The inclusive cross section for i^^ + ^ /i^ + AT is given by 

- — ^ r{xq{x, + 2:(1 - yfq{x, Q^)), (1) 



dxdy TT {Q^+M^y 

where x ~ /2M y = v/ Ey, and v the energy loss in the lab frame, v = E^ — E^. The mass of the nucleon 
is M, and Mw is the mass of the VF-boson, while the Fermi constant, Gp — 1-16 x 10~^ GeV~^. With the 
assumption that the target is an isoscalar nucleon, the quark and antiquark distribution functions, qix, Q^) 
and q{x, Q^), contain valence and sea quark contributions. 

At low energies, the neutrino-nucleon cross section reduces to the usual 4-Fermi approximation and in- 



creases linearly with the neutrino energy. At high energies, ^ M^/2M, the Vt^-propagator in Eq. (1) 
becomes important Q, limiting the growth of Q to (Q^) ^ M'^. For E,j ^ 10^ GeV, the neutrino-nuclcon 
total cross section is dominated by the behavior of the quark densities at a; M^/2MEi, and ^ ^w- 
For neutrino energies above 10^ GeV, the small- a; (i.e., a; < 3 x 10^^) behavior of the quark densities plays 
an important role in determining the cross-section. 

Recently, knowledge of the parton densities at small values of x has been substantially improved by 
the observation of a dramatic increase of the proton structure function, F2{x,Q^) with decreasing x for 
10-"* < a; < 10-2 and 8.5 GeV^ < < 15 QeV^ by the ZEUS and HI Collaborations at HERA g. For 
larger values of , the HERA experiments have measured F2^{x, Q^) down to a value of a; = 2 x 10^^. Thus, 
for neutrino energies above E^ 10^ GeV, the relevant small- a; regime is unconstrained by experiments. 

In the absence of measurements in the appropriate range of (x, Q), the parton distributions are traditionally 
obtained by assuming certain form at some fixed value of Q = Qo ~ few GeV and the distributions are 
evolved to larger values of Q via QGD evolution equations. The parameters in the initial parametrizations 
are determined by fitting a variety of experimental data on structure functions for x > 10^^ and moderate 
values of Q, as well as the recent HERA data Q. 

The theoretical uncertainties in the total cross section for UHE neutrinos in the standard model are 
therefore due to the parameterization at Qo, the extrapolation to values of x below current data and in the 
evolution of the parton distribution functions to Q > Qo- 

Currently there are two theoretical approaches to understanding small- a: evolution with Q, both based on 
perturbative QCD. The first, and more traditional approach is to determine parton densities for Q > Qo by 
solving the next-to-leading order Gribov-Lipatov-Altarelli-Parisi (CLAP) evolution equations numerically 0. 
Crucial to the calculation of the UHE i^N cross section is the small-a parameterization, which in all cases 
has the form for sea quark distributions of ag(a;, Qq) ^ x^^. The second approach to small-a evolution is 
to solve the Balitskii-Fadin-Kuraev-Lipatov (BFKL) equation The solution predicts A ~ 0.5 and a more 
rapid evolution with than the CLAP evolution. 

The parton distribution sets obtained perturbatively using next-to-leading order (NLO) CLAP evolution 
equations are applicable for the calculation of the total neutrino cross section for i?^ < 10^ GeV, and are a 
reasonable starting point for calculations for higher energy neutrinos. The NLO CTEQ3 distributions in 
the deep-inelastic scattering factorization scheme (CTEQ-DIS) have Qo = 1.6 GeV and A = 0.332. These 
distributions are particularly useful because the numerical evolution is provided for a: — > 0, albeit, in a region 
where GLAP evolution is unreliable. 

To estimate the uncertainty in the cross section, we consider two other distribution functions extrapolated 
to small-a. An approximate lower limit is given by the leading order CTEQ3 distribution with the double- 
logarithmic-approximation (DLA) extrapolation below amin ~ 10^*, labeled CTEQ-DLA below. This is an 
approximate solution to the leading order GLAP equations, suitable for A ~ 0. We use a more singular 
distribution than the CTEQ-DIS distributions to indicate the upper range of the cross section. Motivated by 
BFKL dynamics, the MRS D_ set [|l^ is parameterized with A = 0.5. The power A = 0.5 does not change 
significantly with x and Q for small- a; and Q ~ Mw, so we use a power law extrapolation below the limit of 
the numerically evolved distributions, namely for x < Xmin = 10^^. 

In Fig. 1 we present our results for the charged-current cross sections for the three choices of parton 
distribution functions, and compare with the calculation of Ref. which used the Eichten et al. ITT] 
distributions extrapolated using the DLA (EHLQ-DLA). We find our results to be in very good agreement 



Figure 1: The charged current cross section for the CTEQ-DIS, CTEQ-DLA, D_ and EHLQ-DLA parton 
distribution functions. The data point, an average of ZEUS and HI, is from Ref. [12]. 



with the recent measurement of charged-current cross section at HERA IQ . 

The differences in the cross section appear at neutrino energies E^, ~ 10^ GeV, reflecting the different 
behaviors of xqs{x,Q) below x - M^/{2M ■ 10^ GeV)- 3 • IQ-^. The D_ and CTEQ-DIS cross sections 
are essentially equal for ~ 10^ GeV, but the CTEQ-DLA cross section is ^ 15% lower. By E^, — 10^" 
GeV, there is a factor of 2.6 difference between the CTEQ-DLA extrapolation and the D_ calculations of the 
cross section. The CTEQ-DIS result lies a factor of 1.4 above the CTEQ-DLA cross section. At E^, = 10^^ 
GeV, the discrepancy between the CTEQ-DLA and D_ calculations is a factor of 6. One should keep in mind, 
however, that the DLA is an approximation valid for E^y <C 10^^ GeV. 

In order to calculate the number of upward-moving muons that can be detected with future neutrino 
detectors, such as DUMAND II, AMANDA, BAIKAL and NESTOR g, we fold in the muon flux and its 
attenuation due to its passage through the Earth with the probability that a neutrino passing on a detector 
trajectory creates a muon in the rock which traverses the detector. The cross section appears in the exponent 
for the attenuation, and the probability is proportional to the cross section. The details are found in Ref. [|. 
We find that above 10^ GeV, both the muon probability and the neutrino attenuation are very sensitive to 
the small- a; behavior of the sea quark distributions, but their product is not. Thus, for the downward-moving 
muons (and for the contained events) the event rates obtained with MRS D_ or CTEQ distributions are much 
larger than the same with EHLQ structure functions, but the upward rates differ only by 17% for the whole 
energy range. 

We find that future detectors, such as DUMAND II, AMANDA, BAIKAL and NESTOR have a very good 
chance of being able to test different models for neutrino production in the AGNs |^ . For i?^ > 1 TeV, we find 
that 900 upward-moving muons/yr/km'^/sr originating from the diffuse AGN neutrinos, with the atmospheric 
background of 1400 events/yr/km^/sr. For > 10 TeV, signal to background ratio becomes even better, 
the signal of 500 muon events/yr/km^/sr, being almost 20 times the background rate. 
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